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Multiple tracers of the cosmic density field, with different bias, number and luminosity evolu¬ 
tion, can be used to measure the large-scale properties of the Universe. We show how an optimal 
combination of tracers can be used to detect general-relativistic effects in the observed density of 
sources. We forecast for the detectability of these effects, as well as measurements of primordial 
non-Gaussianity and large-scale leasing magnification with current and upcoming large-scale struc¬ 
ture experiments. In particular we quantify the significance of these detections in the short term 
with experiments such as the Dark Energy Survey (DES), and in the long term with the Large 
Synoptic Survey Telescope (LSST) and the Square Kilometre Array (SKA). We review the main 
observational challenges that must be overcome to carry out these measurements. 


I. INTRODUCTION 


Over the next few decades we expect to map out the 
large scale structure of the observable Universe in its 
entirety. By this, we mean that we should be able to 
completely characterize the spatial distribution of mat¬ 
ter from sub-Megaparsec scales out to the cosmological 
horizon, from the present epoch out to the highest ob¬ 
servable redshift, deep into the dark ages of structure 
formation. We will achieve this through a combination of 
redshift (both spectroscopic and photometric) and weak 
lensing surveys allied with intensity mapping and CMB 
experiments. From the spatial distribution of matter, we 
should be able to infer the structure of space-time within 
our cosmic horizon and in doing so dramatically improve 
our understanding of the cosmological model and its fun¬ 
damental assumptions. 

Key to our understanding of the Universe are the very 
largest scales. There we expect to be truly sensitive to 
primordial and general-relativistic aspects of cosmology 
and, one hopes, we should be able to constrain aspects 
of the Universe which are not accessible if we restrict 
ourselves to smaller, Newtonian, scales. The most no¬ 
table aspects, which we will focus on in this paper and 
on which we will elaborate below, are general-relativistic 
effects on number counts and luminosities m, large- 
scale weak lensing and the scale dependence of bias due 
to primordial non-Gaussianity HIT]. 

There have been attempts at constraining the largest 
scales of the Universe with current surveys [5] and of 
quantifying how well we will be able to do so with fu¬ 
ture experiments [Mil. Most of these attempts have fo¬ 
cused on autocorrelations of single tracers, where it has 
become clear that cosmic (or sample) variance is suffi¬ 
ciently large that the resulting constraints are not (or 
will not be) particularly stringent. In other words, the 
fundamentally stochastic nature of the seeds for structure 
combined with the finite volume of any survey severely 
restricts the number of independent modes we can sam¬ 
ple and introduces a large and inescapable floor on the 
uncertainty of any of the fundamental quantities we wish 


to constrain. In m we undertook a detailed, compre¬ 
hensive analysis of the four different types of surveys of 
large scale structure which are currently being planned, 
finding indeed that there were severe limitations to how 
well we could constrain large scale effects. 

An ingenious and alternative approach can be taken if 
one has multiple tracers (such as galaxies or line-emitting 
species) with different bias and luminosity functions. As 
was first shown in |lfij and then further developed in 
m, with multiple tracers it is possible to effectively ” di¬ 
vide out” the stochastic dependence of the cosmological 
density field allowing us to beat cosmic variance for a 
few, specific, observables such as the growth rate and 
clustering bias. While this approach was originally ad¬ 
vocated to target the scale dependent bias of primordial 
non-Gaussianity m, we expect that such an effect nat¬ 
urally extends to the general-relativistic corrections in 
large-scale structure surveys. 

In this paper, we further flesh out the idea of using 
multi-tracer methods to constrain what we dub ultra 
large-scale effects. We focus on surveys that have two 
key properties. First, they must be able to cover large 
volumes of space with an appreciable number density of 
sources. Second, we need to be able to pick out two or 
more different tracers with a sufficiently high bias but 
also with sufficiently different evolution properties that 
the multi-tracer technique can be effective. In subsec¬ 
tion III B| we derive a simple expression which can be 
used as a guide for what the best combinations of tracers 
(and their properties) are to efficiently constrain general 
relativistic effects. Both of these requirements lead us 
to choose photometric redshift surveys (focusing on DES 
and LSST), radio-continuum surveys and intensity map¬ 
ping surveys (primarily from the SKA). Indeed we find 
that a combination of these surveys will give us remark¬ 
able constraints on ultra large-scale effects. 

We structure our paper as follows. In Section jlT] we re¬ 
mind ourselves of the various effects that come into play 
on large scales, explaining their origin and how they af¬ 
fect the observables we will consider. In Section IIIII we 
describe the Fisher forecasting formalism we will be using 
and employ it to derive a simple expression which shows 
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how the multi-tracer method works to our advantage, for 
isolating the general-relativistic corrections to large-scale 
surveys. This expression gives us a few pointers to what 
the optimal surveys will be for us to achieve our goals. 
In Section IV we focus on optical and near infra-red pho¬ 
tometric redshift surveys, specifically on LSST, in Sec¬ 
tion |V] we look at cosmological radio surveys, specifically 
continuum and intensity mapping as will be carried out 
by the SKA, and in Section VI we will see how cross- 
correlating these different surveys will contribute to our 
overall goals. Our results are truly promising and in Sec¬ 
tion m their limitations and future prospects are pre¬ 
sented. In particular we speculate on how precise mea¬ 
surements on ultra-large scales will contribute to our un¬ 
derstanding of the fundamental cosmological model as 
well as how they might help us further characterize the 
large scale properties of the Universe. 


II. COSMOLOGICAL OBSERVABLES ON 
ULTRA-LARGE SCALES 

The aim of this section is to introduce the two main 
observables on ultra-large scales that we will study here, 
as well as to establish the notation that will be used in 
the rest of the paper, which follows the one used in m- 
More thorough discussions of these topics can be found 
in, e.g. 0 for relativistic effects and [7] for primordial 
non-Gaussianity. 


A. Relativistic effects in large-scale structure 


there are other non-zero contributions at linear order that 
potentially contain information about the nature of the 
gravitational interactions perturbing the photon trajec¬ 
tories. The complete set of terms at linear order was 
presented by [S] and H], and can be classified into four 
main contributions: 


AAr(z, n) = 


d\nN' 

-j—Sr,+Sii+2S± 


d\nN' 
9 In L* 


( 1 ) 


where A 7 (? 7 ,lnL*) is the physical number density of 
sources at conformal time 77 with luminosities above L* 
(i.e. the cumulative luminosity function), related to the 
survey’s flux limit Fcut via 

L*(z) = 47r(l-h Fcut- (2) 

The four terms on the right hand side of Eq. corre¬ 
spond to the intrinsic perturbation in the number den¬ 
sity of sources {Sj^f), the difference between the source’s 
conformal time 77 and the conformal time inferred from 
the redshift 77 * (Sr]), and the corresponding differences in 
the radial distance ((5||) and angular diameter distance 
(i5_l) of the sources. As discussed above, RSDs and lens- 
ing magnification contribute to Jy and (5 _l respectively. 
In a given gauge, we can expand each of these in terms 
of the metric, density and velocity perturbations. This 
has been done explicitly in the literature, and we present 
expressions for the complete set of terms in Appendix [A| 
Note that the amplitude of the observed perturbation 
in the number counts as a function of redshift depends 
on three tracer-specific quantities: 


The aim of large-scale structure experiments is to map 
out the three-dimensional distribution of the matter den¬ 
sity field in the Universe. This is most commonly done 
using the number counts of astrophysical sources in the 
sky, using their angular positions and redshifts as prox¬ 
ies for their three-dimensional coordinates. This means 
that we use the redshift and direction of the photons 
emitted by these sources in order to infer the comoving 
volume of the patch over which they are distributed; any 
perturbation in the trajectory or the frequency of these 
photons along the lightcone, will affect our volume esti¬ 
mates, and consequently the measured three-dimensional 
density map. 

Two perfect examples of such effects are redshift- 
space distortions (RSDs), by which the emitted photons 
are further redshifted by the peculiar velocities of the 
sources, and gravitational leasing magnification, the per¬ 
turbation in the trajectory of the photons caused by the 
gravitational field of the intervening matter. RSDs dis¬ 
place the sources in redshift-space coherently along the 
line of sight, while leasing magnification modifies the 
solid angle subtended by a cluster of sources. While these 
two effects correspond to the two largest contributions to 
the total fluctuation in the observed number counts, after 
the intrinsic perturbation in the source number density. 


• The clustering bias 6 , relating the intrinsic per¬ 
turbation in the number density with the fluctu¬ 
ations in the comoving-gauge matter perturbation, 
defined in Fourier space as [SmillS]: 

(3) 

where is the matter perturbation in the syn¬ 

chronous comoving gauge and v is the peculiar ve¬ 
locity in the Newtonian gauge. 


• The slope of the background number density of 
sources as a function of time, which is now cus¬ 
tomary to encode in the so-called evolution bias: 


_ d\n{o? M) 
91n(a) 


(4) 


• The slope of the number density of sources as a 
function of intrinsic luminosity, usually written in 
term of the magnification bias: 


91n(A7) 
2 d\nL 
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In general, the astrophysical observable containing the 
most cosmological information is the perturbation in the 
number counts in a set of redshift bins 

^i(n) = J dzWi{z)A{z,h), (6) 

where Wi{z) is the the window function defining the i—th 
redshift bin (normalized to 1 when integrated over the 
whole redshift range). To make matters simpler, we will 
present our forecasts in terms of the harmonic coefficients 
of Ai(n), defined as 

Ofm = J dhA,{h)Ye^{h). (7) 

The main advantage of using the a^mS is that statistical 

isotropy makes them uncorrelated, so that: 

^rnm' ■ ( 8 ) 

This equation also defines the angular power spectrum 
containing all of the possible cross-correlations be¬ 
tween two redshift bins i and J. 

As shown in m, the angular power spectrum can be 
computed as 

7°° lik 

-^V{k)Al{k)Al{k), (9) 


to the first three. Moreover, since accurate clustering 
measurements on these scales can be severely hampered 
by observational systematics, as we will review later, de¬ 
tecting the effect of these terms is very challenging. In 
particular, as was shown in m, the amplitude of these 
terms is too small to be observable with any significance 
using any of the main cosmological probes targeted by 
next-generation large-scale structure experiments indi¬ 
vidually. The main reason for this is the large cosmic 
variance on these scales, both due to the small number 
of super-horizon modes available at redshifts 2 : < 3 and 
to the larger amplitude of matter fluctuations at late 
times. In order to detect these relativistic corrections, 
it is therefore necessary to reduce the statistical noise 
below the cosmic variance limit, a feat that might be 
achievable using the multi-tracer techniques that we re¬ 
view in section III A We will quantify the detectability 
of these GR corrections as was done in m. by defining 
an effective parameter cgr multiplying the terms 

and We do not include the lensing mag¬ 

nification term A^ under the umbrella of cgr, since this 
effect has already been detected observationally [Ml El] , 
and its contribution is also non-negligible on small an¬ 
gular scales. However, for completenes, and given that 
gravitational lensing is clearly a General-Relativistic ef¬ 
fect, we will also forecast for its detectability by defining 
a separate amplitude parameter for it, e-wL- 


where V{k) is the dimensionless primordial power spec¬ 
trum, which we assume to be close to scale invariant: 
V{k) = As(fc/fco)""’“^, and A\ is a transfer function con¬ 
taining the terms in Eq. [l] projected along the radial 
direction in each redshift bin. In the Newtonian gauge 
A\ can be expanded as a sum of 10 individual terms, 

A = A° -f -I- A'" -I- A"^^ -I- A^^ 

-f API-h AP2 + (10) 

where the complete expressions for each term are pre¬ 
sented in Appendix fAj 

The first three terms, A^, AP^p* and A^ are the dom¬ 
inant contributions to the total fluctuation, and corre¬ 
spond to the intrinsic perturbation in the source number 
density, the Kaiser-term for redshift-space distortions 
and the effect of lensing magnification respectively. The 
remaining 7 terms are subdominant and have never been 
accounted for, nor detected in the analysis of any dataset. 
Of these, the terms A^^ and A^^ are extra RSD effects 
due to evaluating the background terms at the redshift 
perturbed by the Doppler effect, and the remaining terms 
correspond to similar effects caused by gravitational red- 
shifting and lensing (e.g. the Shapiro time delay AP^ 
or the integrated Sachs-Wolfe effect |22l ESj). A 

great deal of the work presented in this paper is devoted 
to studying the detectability of these last 7 terms, and 
the potential benefits that such a detection could entail. 

The contribution of these terms to the observed fluc¬ 
tuations is most relevant on the scale of the cosmic hori¬ 
zon, and even there they are subdominant with respect 


To end this section, it is worth noting that, even 
though the form of the fluctuation Ai(n) has been mo¬ 
tivated for experiments targeting a flux-limited popula¬ 
tion of discrete sources, it is straightforward to extend it 
to the case of line-emission tomography (most com¬ 
monly called intensity mapping ). The main difference 
between these two cases is that, while in a survey of dis¬ 
crete sources the proxy for the density perturbation is 
the fluctuation in the number of objects detected above 
a given flux in different patches of equal size, the ob¬ 
servable in intensity mapping is the combined emission 
from all of the sources in these patches. Since lensing 
effects perturb angular separations in the same way as 
source luminosities, the net effect is an exact cancella¬ 
tion of the linear-order perturbation in the angular di¬ 
ameter distance for intensity mapping (i.e. the i5_l term 
in Equation]^. Thus, for the purposes of predicting the 
cosmological signal, intensity mapping experiments can 
be treated as a survey of discrete sources with a magni¬ 
fication bias given by the critical value sim = 2/5, and 
with an evolution bias given by: 


IM ^ d\npa 

aina ’ 


where Pa is the background comoving density of the line- 
emitting species under study. 
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B. Primordial non-Gaussianity 


One of the fundamental assumptions of the simplest 
versions of the standard cosmological model is an almost- 
negligible amount of non-Gaussianity in the statistics of 
the primordial density field m- A simple (although not 
universal) formalism [ 28 ] to quantify deviations with re¬ 
spect to this premise is to decompose the primordial grav¬ 
itational potential into a linear and a quadratic func¬ 
tion of a Gaussian random field (j)G- 

^ = (j)G + /nl {(I^g - {(fc))- ( 12 ) 


Since the matter density field gradually becomes non- 
Gaussian due to the non-linear nature of gravitational 
collapse, the level of primordial non-Gaussianity is eas¬ 
ier to determine using early-Universe probes, such as the 
GMB, which are well inside the linear regime. Gurrent 
constraints based on higher-order GMB statistics favour 
a rather small level of non-Gaussianity, with |/nl| 7 

However, other statistical techniques have been de¬ 
veloped with the aim of improving these constraints us¬ 
ing low-redshift three-dimensional datasets. In particu¬ 
lar it has been noted that a non-zero primordial non- 
Gaussianity would affect the formation and therefore the 
clustering statistics of biased tracers of the matter distri¬ 
bution laiz]. The effect can be encoded as a modification 
of the Gaussian clustering bias bciz) —>• &g(z)-|-A6 (z, fc), 
with a scale dependence given by 


A6(2, k) = 3/nl 


[bciz) - Ij^M-gg^e 
k^T{k)D{z) 


( 13 ) 


corrections in terms of egr and for the level of primor¬ 
dial non-Gaussianity parametrized by /nl- We do so us¬ 
ing the Fisher matrix approach, extending the formalism 
that was used in m to include different tracers. 


A. Fisher matrix with multiple tracers 


Let us consider an experiment targeting a number of 
tracers, each labelled by a Greek index a; each tracer 
results in a number of maps of the sky, labelled by a Ro¬ 
man index i. For the purposes of this paper, these maps 
correspond to different redshift bins or frequency bands. 
All the cosmological information is contained in the har¬ 
monic coefficients for each tracer and map, which we label 
by a combined index (a, i) (i.e. are the harmonic 

coefficients for the a-th tracer in its i-th redshift bin). 
For Gaussian random fields, the most important observ¬ 
able is the power spectrum, defined as the covariance of 
the 

(is) 

We will assume that all tracers can be cleanly sepa¬ 
rated into signal and noise contributions -I- 

both components are statistically inde¬ 
pendent. Furthermore, we will assume that the noise is 
uncorrelated between different tracers: 

( 16 ) 


where Sc — 1.686 is the linearized collapse threshold, 
T(k) is the matter transfer function and D{z) is the linear 
growth factor. 

This scale-dependent bias induces extra power on large 
scales with a ~ scale dependence, and it is easy to 
show that this extra power should be more evident on 
scales or the order of the horizon /cng /nl-TIo- Thus 
this signature of primordial non-Gaussianity presents 
characteristics that are very similar to the signal of the 
GR corrections presented in the previous sections (in fact 
it has been suggested that they could potentially be de¬ 
generate [H [201 E]), and therefore it makes sense to 
study both effects simultaneously. 

Including our three main observables (/nL) cgr and 
ewl), the total number count fluctuation is: 

Af = Af (/nl) + Af®° -I- ewL A^ -I- cgr Af (14) 

where every term is implicitly scale and time-dependent, 
and A^^ = A^^ -h A^^ -f A^^ -h A^^ -H A^^ -1- A^"^ -H A^®'^. 

III. FORECASTING FORMALISM 

The aim of this paper is to produce forecasts simul¬ 
taneously for the detectability of the general relativistic 


This assumption would not be correct, for instance, if 
the two tracers were two different but not disjoint galaxy 
populations. For most applications it is also safe to as¬ 
sume that the noise power spectrum is uncorrelated be¬ 
tween different redshift bins. 

Arranging the combined indices (a, i) as a single one¬ 
dimensional index, we can write the and the power 

spectra as a vector and a matrix respectively. For 
example, for two tracers a = {1,2} distributed over 2 
and 1 redshift bins respectively, we would have a-im = 

^ power spectrum 




o o 

to to 

V 

> / 


where we have omitted the redundant components 
strictly below the diagonal. It is then easy to show that 
the Fisher matrix for the observable is the same as for 
the single tracer case, where the power spectrum matrix 
now contains also all possible cross-correlations between 
different tracers. Thus, for a set of parameters { 0 ;^}, the 
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covariance matrix = ((0^ — 0^){0i, — 9^)) can be ap¬ 
proximated by the inverse of where 

o/ -I- 1 

F^u = /sky 51 ^^Tr [a^Q C.-ia.Q C7I] ( 18 ) 

e=2 


and is the partial derivative with respect to the /r-th 
parameter. In our forecasts the derivatives were com¬ 
puted using central finite differences, 




f{9^ + Se^) - f{9^ - 59^) 
269 ^ 


0{69^), 


(19) 


where the intervals 69^^ are such that the derivatives con¬ 
verge to the required numerical accuracy. 

As we pointed out above, we can divide the power 
spectrum up between two uncorrelated signal and noise 
contributions, -|- N^. The cosmological power 

spectra Cf are given by Equation and were calcu¬ 
lated using a modified version of the public CLASS code 
[iniiii]- This version is based on the code used in HU, 
with added capability to handle an arbitrary number of 
tracers m- For our fiducial cosmology we adopted a 
model consistent with the best-fit flat ACDM parame¬ 
ters from Planck [S3], given by (JIm, fb,h,w,As,ns) = 
(0.315,0.156,0.67,-1,2.46 x 10-9,0.96). We further set 
/nl = Oj the value for Gaussian initial conditions, and 
ecR = ewL = 1- In our forecasts we use these 9 free 
parameters, as well as a number of nuisance parameters 
associated with the bias functions of the tracers under 
study, which we describe below. Any priors on these pa¬ 
rameters, correlated or not, were added as Gaussian pri¬ 
ors to the bare Fisher matrix (F^^ —>• + (Gp™'')”/). 

In particular for the cosmological parameters listed above 
(excluding /nl, ecR and cwl) we used the covariance 
matrix estimated from the corresponding Planck MCMC 
chains. We must also note that, since cgr and ewL are 
only effective parameters, used to characterize the level 
of detection of the GR and lensing effects, when forecast¬ 
ing for /nl we fix both cgr and cwl to their actual value 
of I instead of marginalizing over them. 

The uncertainty with which a given parameter can be 
constrained usually depends on the largest and smallest 
accessible scales. The smallest angular scale included in 
our forecasts is given by in Eq. [T^ while the smallest 
radial scale is given by the width of the redshift bins used 
in the analysis. These scales should be chosen on the ba¬ 
sis of the angular and radial resolution of the experiment 


Tracer 

bins 

Bin edges 

Int. map. (SKAl) 

4 

0.0, 0.7, 1.3, 2.0, 3.0 

Cont. survey (SKAl) 

4 

0.0, 0.5, 1.0, 2.5, 5.0 

Photo. (LSST), red 

4 

0.0, 0.35, 0.7, 1.05, 1.4 

Photo. (LSST), blue 

5 

0.0, 0.5, 1.0, 1.5, 2.2, 3.0 

Photo. (DES), red 

4 

0.0, 0.3, 0.6, 0.9, 1.2 

Photo. (DES), blue 

4 

0.0, 0.4, 0.8, 1.2, 1.5 


TABLE I: Redshift bins used for the bias function nuisance 
parameters, for each tracer. 



0 1 2 3 4 5 


z 

FIG. 1: Clustering bias (upper panel), magnification bias 
(middle panel) and evolution bias (lower panel) for the dif¬ 
ferent tracers considered in this work: intensity mapping (or¬ 
ange), first continuum sample, made up of star-forming galax¬ 
ies, starbursts and radio-quiet quasars (cyan), second contin¬ 
uum sample, made up of FRI-type AGN (green), combined 
continuum sample (grey), red LSST galaxies (red) and blue 
LSST galaxies (blue). In each curve, the markers show the 
nodes used as free parameters in the marginalization over 
these bias functions. Notice the very different evolution biases 
for red and blue galaxies. 


(e.g. the beam size for intensity mapping or the red¬ 
shift uncertainty for photometric surveys), or in terms of 
the smallest scales for which we can trust our theoretical 
prediction (e.g. the non-linear scale). Fortunately, since 
the observables we are interested in are mostly relevant 
on ultra-large scales, our forecasts are relatively insensi¬ 
tive to the choice of minimum scales. By default we use 
^max = 500 for all the experiments considered here, and 
address the redshift binning scheme on an individual ba¬ 
sis. We further assume that all these experiments will be 
able to probe angular scales down to £min = 2, and we 
asess how the constraints degrade as a function of imin- 
We treat the uncertainty in the bias functions 
{b{z), s(z) and fevoiz)) using the same method that was 
used in m- the mean value for each function is com¬ 
puted in a small set of wide redshift bins, and we use 
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the linear interpolation between these values as the fidu¬ 
cial bias function passed to CLASS. We then treat each of 
the aforementioned mean values as free parameters that 
we include in our computation of the Fisher matrix and 
marginalize over. The redshift bins used to define these 
bias parameters for the different probes considered here 
are given in Table |Tj and the resulting bias functions are 
shown as a function of redshift in Figure 

B. Why multi-tracer surveys beat cosmic variance 

It was first noted in m that using different tracers 
deterministically related to the same density field it is 
possible to measure certain parameters with uncertainties 
smaller than the cosmic-variance limit. The aim of this 
section is to briefly describe the mechanism by which this 
is possible, and how we could benefit from it in order 
to detect the relativistic corrections presented in Section 

ins 

Let us consider a simplified scenario, where we mea¬ 
sure the fluctuations in the source number counts of two 


different tracers (a = 1, 2) in the same single redshift bin. 
For the sake of the argument let us assume that only two 
terms contribute to the total perturbation: 

aim = ba AL + U e ^fm + ■ (20) 

Here A'^ and A*^ are the terms corresponding the intrin¬ 
sic perturbation in the matter density and a generic rel¬ 
ativistic correction respectively, ba is the clustering bias, 
fa is a tracer-dependent bias function affecting the GR 
correction (i.e. similar to the evolution or magnification 
biases) and is a noise term. We have explicitly ex¬ 
tracted the amplitude of the GR term as a parameter, e, 
equivalent to ecR- This can be rewritten as 

a7m = {be t! + fa e Tf) (21) 

where Tf and are the transfer functions correspond¬ 
ing to the density and GR correction terms, and is 
the primordial perturbation. 

The covariance matrix for both tracers in each multi¬ 
pole is 


C = 


/ hlV{T^ + 

\b\V{T^ + - 


-Ni 

/32eTG) 




blVijT^ + P2eT^)^ 


P2eT^) \ 

bAa ) 


( 22 ) 


r 


where V is the power spectrum of 7 = fo 2 /^ii Pi = 
fi/bi, p 2 = f 2 /bi, Na is the noise power spectrum for 
each tracer and we have assumed that ( 711712 ) = 0. Note 
that we have omitted the dependence on the multipole 
order i. 

We can estimate the best-case uncertainty on the am¬ 
plitude e by computing the inverse of its Fisher matrix 
element a‘^{e) ~ l/Feg. Using = ^Tr [(C“^i9gC)^] we 
obtain 


R /3i6 

2 

-U 

'yi? + P 2 I 

_p2 - lPl_ 

-f yvi 

. P 2 - iPi_ 


where we have defined the ratios Xa = Na/iV {biT^Y) 
and R = /T^ ^ and we have linearized the result with 

respect to Several conclusions can be drawn from 
Eq. 

1. The error on e is completely limited by the noise 
level, and not by cosmic variance. This could 
have been anticipated by looking at Equation 

in the noiseless case (71“ —)■ 0): by taking the ratio 
a\mlaim it would be possible to cancel the stochas- 
ticity of <1)° and determine ba , fa and e free of cos¬ 
mic variance. This is the key to understanding the 
multitracer technique. 

2. If the amplitude of the relativistic correction is sub¬ 
dominant, (i.e. R ^ 1), as is the case for the terms 


studied here, the uncertainty in e becomes propor¬ 
tional to i?, and thus achieving a{e) 1 would 
still require very low noise levels and/or appropri¬ 
ate values for the bias parameters fa- 


3. The effect of the factor P 2 — iPi in the denomi¬ 
nator is twofold. Firstly, by choosing two tracers 
with very different bias functions (e.g. /2 S> /i) it 
is possible to further reduce the uncertainty on e. 
This is key to understand the results presented in 
subsequent sections. Secondly, in the case Pi = P 2 
and 7 = 1 the uncertainty on e diverges, which is a 
consequence of the fact that in this case both trac¬ 
ers are actually the same one, and nothing can be 
gained from treating them as being different. 


Although this result is based on a toy model, it allows 
us to build a picture regarding the optimal tracers that 
should be used in order to achieve a measurement of the 
GR corrections: the ideal case would be to combine mea¬ 
surements of two (or more) high signal-to-noise tracers 
with very different magnification and/or evolution biases. 
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IV. OPTICAL AND NEAR-INFRARED 
EXPERIMENTS 


A. Photometric redshift surveys 

The observational study of cosmological large-scale 
structure has so far been dominated by galaxy surveys 
performed in optical and infrared wavelengths |34H36j . 
In the ideal case, redshifts and angular positions are 
measured for a homogeneous sample of targets, and the 
three-dimensional matter density distribution is probed 
through the fluctuations in the galaxy number density. 
Accurate redshift measurements are, however, very time 
consuming, especially for faint sources, and therefore 
spectroscopic redshifts are usually only available for a 
subsample of all the photometrically-detected sources. 
This has two critical drawbacks: firstly, it prevents spec¬ 
troscopic redshift surveys from covering the largest ac¬ 
cessible scales, which is key to studying the observables 
described in section |IT] Secondly, the sub-sampling of the 
galaxy population produces relatively high shot-noise lev¬ 
els, which dominate the error budget for parameters like 
/nl and cgr in the multi-tracer scheme. These prob¬ 
lems can be partially circumvented using photometric 
redshifts (photo-zs), i.e. redshifts inferred from the ob¬ 
served source fluxes in a number of broad bands. On the 
one hand, these redshift estimates are much less accurate 
than their spectroscopic counterparts, which limits the 
ability of photometric redshift surveys to study cluster¬ 
ing statistics on small radial scales. On the other hand, 
photo-zs can, in principle, be estimated for all the sources 
detected above a given flux threshold, and thus the prob¬ 
lems of volume coverage and shot noise are greatly allevi¬ 
ated. Since the observables studied here are only relevant 
on the largest scales, and their measurement can be op¬ 
timized by reducing the noise levels of the experiment, 
photometric redshift surveys seem to be ideally suited 
for our purposes. 

We will focus our discussion on the Large Synoptic Sur¬ 
vey Telescope (LSST, [37]), since it represents the widest 
20000deg^) and deepest (r ~ 27) photometric survey 
planned in the foreseeable future, but we will also present 
forecasts for the ongoing Dark Energy Survey |38) . 

In this section we focus on the advantages of combin¬ 
ing clustering information for different galaxy samples 
observed within the same experiment. In particular we 
will use a simplified picture in which we assume that the 
full galaxy population can be separated into two disjoint 
samples of “red” and “blue” galaxies: 


pears only after the peak in the star formation 
rate (z ~ 1-9), and its comoving number density is 
known to decrease sharply at redshifts larger than 
~ I. 

The simplest models of galaxy formation suggest 
that the earliest galaxies should have formed in the 
peaks of the density held. Thus the red popula¬ 
tion can be associated with larger-mass halos and 
a larger clustering bias. 


• Blue galaxies, on the other hand, are character¬ 
ized by a population of young stars, showing strong 
Balmer emission lines. They are usually associated 
with a spiral morphology, smaller halo mass and 
lower clustering bias. The number density of blue 
galaxies shows no sharp increase or decrease, and 
LSST should be able to observe them up to much 
higher redshifts than the red population. 


This simplihed model of the galaxy population is based 
on the observed bimodal nature of the color-magnitude 
diagram [3^. The actual picture is nowhere as simple 
(e.g. red spirals and blue ellipticals are known to exist), 
but this model allows us to study the potential improve¬ 
ment in the figure of merit for ultra large-scale observ¬ 
ables caused by combining clustering observations for two 
populations with similar characteristics. 

Of particular interest would be the sharp decay in the 
number density of red galaxies. As was shown in m, 
this is associated with a large evolution bias, and en¬ 
hances the amplitude of the relativistic corrections de¬ 
scribed in section IIA although this enhancement was 
shown to be too small to grant a detection of these ef¬ 
fects using single-tracer observations. However, as dis¬ 
cussed in section Enl combining red and blue galaxies, 
the latter having a much lower evolution bias, could po¬ 
tentially boost the figure of merit for cgr- 

As was done in m, we base our estimates for the num¬ 
ber densities of both populations, as well as their mag¬ 
nification and evolution biases on luminosity functions 
estimated by SOI mi in the r-band. The details of this 
model are summarized in |15j . Regarding the clustering 
bias, we parametrize its redshift dependence for the red 
and full galaxy samples as: 


brcd{z) = l + z, bfuiiiz) = 1 + 0.84 z, (24) 

based on the simulations of [42] and the measurements of 
[33] . We then assume that the clustering bias of the full 
sample is a weighted average of the red and blue samples, 
and compute the clustering bias of the latter as: 


• Red galaxies are usually identihed with an early 
population of galaxies where star formation has 
ended. This produces a red spectrum, usually char¬ 
acterized by a deep Balmer/4000A break. This 
spectral type is also commonly associated with 
an elliptical morphology, supported by random 
motions rather than angular momentum, possi¬ 
bly caused by mergers. The red population ap- 


, , , n{ull{z)biun{z) - nred{z)bred{z) 

^hluey^j — _ / \ • 

Finally we assume a magnitude cut of f = 25.3 (r ~ 
26), corresponding to the so-called “gold sample”, and 
a sky fraction /sky = 0.5 20000 deg^) for LSST. We 

find that the total number counts as well as the mag¬ 
nitude distribution for each galaxy type obtained using 







FIG. 2: Redshift distribution (thick solid black lines) and radial window functions (thin coloured lines) for the redshift bins 
used here for the red (left panel) and blue (right panel) galaxy populations for LSST. 


these models is in good agreement with previous esti¬ 
mates [ 371111 ]. For the Dark Energy Survey we will as¬ 
sume a magnitude limit r = 24 and a sky fraction of 
0.125 (- 5000deg^). 

The redshift binning scheme was determined for both 
samples in terms of their expected photometric redshift 
errors. We assume Gaussianly distributed errors with a 
variance crz(z) = uq{1 + z). The photo-z requirements for 
LSST as quoted in [ 37 ] are (Tq < 0.05 with a goal of 0.02. 
Since the spectral characteristics of red galaxies, partic¬ 
ularly the prominent Balmer break, make their photo-zs 
more reliable, we assumed that the goal quoted above 
will be attained for this sample, and we conservatively 
set the photo-z error parameter for blue galaxies to the 
requirement limit. I.e.: 

= 0.02, crj?'"® = 0.05. (26) 


errors assumed here, w takes the form 


w{z\zo,Zf) = 


erf ( - erf f ^ 


■ (28) 


Using this prescription we divide the red and blue sam¬ 
ples into 15 and 9 redshift bins for LSST (10 and 7 bins 
respectively for DES), shown in figure together with 
their overall redshift distribution. 

The noise power-spectrum is given by = Sijjn'^, 
where 


dzN{z)w\z) 


(29) 


is the angular number density of sources in the i-th red¬ 
shift bin with a photo-z window function in*, given by 
Eq. 


For DES we kept = 0.05 but increased the error 

parameter for red galaxies to ctq®'^ = 0.03, since the qual¬ 
ity of the photo-zs should be degraded by the absence of 
the u-band in the DES filter set. We define the redshift 
bins used for each sample as top-hat bins in the space of 
photo-zs with a width three times the size of the photo-z 
uncertainty at the median redshift of the bin. This width 
was chosen in order to reduce the correlations between 
redshift bins caused by the tails of the photo-z proba¬ 
bility distribution. Thus, a top-hat bin in photo-z space 
with edges zq and z/ has a true-z window function given 
by 

lU(z) oc W(z)'u;(z|zo,z/), (27) 

where w{z\zo, Zf) is the integral of the photo-z distri¬ 
bution over the redshift bin. Eor the Gaussian photo-z 


B. Forecasts 

We begin our discussion with an assessment of the de¬ 
tectability of the GR effects in the simplest and most 
optimistic case: we marginalize only over cgr and /nl 
(since the latter could potentially be very degenerate with 
the former) and fix all other parameters to their fiducial 
values. This would represent what an experiment aiming 
at a first detection of the GR corrections would do. The 
results using the red sample only, the blue sample only 
and the combination of both (in a multi-tracer sense) are: 

O'red(eGR) = 1.531, Crblue(eGR) = 3.094, 
f^joint{cGR) ~ 0.095. 

We can see that the combination of both samples causes 
a dramatic improvement in the figure of merit. Erom our 
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FIG. 3: Improvement of the detection level for the relativis¬ 
tic corrections as a function of Fred (see Eq. 301, which 
parametrizes the amplitude of the difference in the values of 
/evo and s between red and blue galaxies. As discussed in 
Section mE the uncertainty on egr decreases as the values 
of these two parameters for both populations become more 
different. 


discussion in section [ill B[ we can a priori attribute this 
to two causes: the low shot-noise levels for LSST and the 
very different bias parameters of blue and red galaxies, 
especially the evolution bias, as can be seen in Fig. 

In order to show this explicitly we have carried out the 
following exercise. We recomputed our forecasts for a 
fictitious LSST experiment with an ordinary blue sample 
and a red sample with the same number density and clus¬ 
tering bias shown in Figures and [l] but with magnifica¬ 
tion and evolution biases given in terms of a parameter 
Tied Q-S: 


^ + (1 - F;ed)s““^ 

fefo ^ T,ed/e^:o" + (1 ” T.ed)/e“o^ (30) 

For a value of Fled = 0 both samples would have the 
same magnification and evolution biases and, according 
to our discussion in Section jlUBj the uncertainty on cgr 
would not suffer any significant enhancement. As we in¬ 
crease the value of F^ed towards 1 we should gradually 
recover the result quoted above. This is explicitly shown 
in Figure]^ which presents the forecasted uncertainty on 
cgr as a function of F^ed- A 2cr detection of egr would 
already be achievable, in this optimistic setup, just with 
Fred 0.2. 

Since it seems like it could be possible to detect the 
relativistic corrections using the multi-tracer method, we 
must now explore the possible degereracies that could, 
on the one hand, preclude this detection and, on the 
other, affect the measurements of other parameters. We 
do not observe strong degeneracies between egr and any 


cosmological parameter, and only a very mild one with 
/nl- The situation is different, however, for the nuisance 
bias parameters. Since the amplitude of the relativistic 
corrections is partially determined by the values of the 
magnification and evolution biases of the samples under 
study, it is probably very difficult to claim a detection 
of these effects without good prior information on s and 
/evo- We have studied the dependence of the uncertain¬ 
ties on /nl, ecR and ewl on a Gaussian prior, constant 
in redshift, for these two bias parameters. The results 
are shown in Fig. j^in the same three cases explored be¬ 
fore (red and blue samples independently and jointly). As 
should be obvious, the detection level of the lensing mag¬ 
nification term depends critically on the precision with 
which the magnification bias is known, and an error on 
s{z) of the order of ^ 0.01 would be necessary to opti¬ 
mize this measurement. On the other hand both /nl and 
EGR are almost insensitive to the prior on s. The situa¬ 
tion is different for fevo- both /nl and egr are sensitive 
to the prior on fevo, although only a very broad prior of 
A/evo = 0(1) would be necessary to optimize their mea¬ 
surement in the single-tracer case. For the multi-tracer 
analysis, however, the detection of the GR corrections de¬ 
pends critically on having precise enough measurements 
of the evolution bias, which as we saw is responsible for 
significantly boosting the signal. The optimal case would 
correspond to a precision of A/evo ^0.1. We have also 
studied the possible dependence of the uncertainties on 
our key observables with the priors on the clustering bias, 
finding no significance degeneracies. 

The availability of such tight priors on the evolution 
and magnification biases will depend directly on our abil¬ 
ity to model the joint redshift and magnitude distribution 
of the galaxy samples under study. This could be chal¬ 
lenging for photometric redshift surveys, given the poor 
redshift uncertainties, and would probably entail relying 
on external datasets or on a spectroscopic follow-up of a 
representative subsample of the survey. In order to gain 
some intuition regarding the feasibility of these measure¬ 
ments we have estimated the uncertainty on the evolution 
bias of red galaxies (responsible for the enhanced ampli¬ 
tude of the GR effects) given the uncertainties in the 
measurements of the luminosity function of red galax¬ 
ies of [40] on which our model is based. In the high- 
redshift end of the distribution we obtain an uncertainty 
A/evo 2, which decreases towards smaller redshifts as 
the number counts increase. In the future these uncer¬ 
tainties should improve significantly due to the availabil¬ 
ity of larger and deeper datasets, and therefore we will 
present our forecasts in a pessimistic and an optimistic 
scenario, where we will assume the priors: 

Pessimistic —J-As = A/evo = 1, 

Optimistic —J-As = A/evo = 0.1. 

Figure [^presents the joint forecasted uncertainties on 
/nl , CGR and ewl for LSST in the optimistic scenario for 
the three cases mentioned above: red-only (red), blue- 
only (blue) and joint analysis (black). The numerical 
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FIG. 4: Dependence of the uncertainty on /nl, egr and ewL on the assumed Gaussian constant prior on the magnification bias 
(left panel) and the evolution bias (right panel). The significance of the detection of the lensing magnification term depends 
very strongly on the prior on s. Likewise the error on cgr is particularly sensitive to the prior on /evo when both tracers are 
combined, and the detection of GR effects only becomes optimal for A/evo < 0.1. 



FIG. 5: Joint constraints on /nl, egr and ewL for the red 
and blue samples (red and blue contours respectively) as well 
as for a joint analysis of both (black contours) for optimistic 
priors A/evo = As = 0.1. 


values for these uncertainties are given in Table [IT| Note 
that, unlike in the case of ecR there is only a mild im¬ 
provement in the figure of merit for /nl and ewL in the 
multi-tracer analysis with respect to the deepest LSST 
blue sample alone. The main reason for this is that, un¬ 
like in the case of the evolution bias, the differences in 


the magnification and clustering biases of both tracers 
are not so large. 

We have also produced forecasts for the Dark Energy 
Survey, using the same models adopted for LSST with 
a magnitude limit of r = 24 and /sky = 1/8. The re¬ 
sults are also included in Table [IT} DES should, in the 
best-case scenario, be able to make a ~ Scr detection of 
the relativistic corrections. This result, however, could 
be compromised by the possible systematic effects that 
could dominate the clustering statistics on large angular 
scales. We will discuss these in Section Eul 


V. RADIO EXPERIMENTS 
A. Cosmological radio surveys 

With the forthcoming wide-area radioastronomy facil¬ 
ities, the field of observational large-scale structure will 
soon begin to reap the benefits of observing in the ra¬ 
dio regime. The low atmospheric absorption and dust 
obscuration in a wide range of radio frequencies makes 
it possible to observe objects at significantly higher red- 
shifts than are usually targeted in optical/NIR surveys, 
and in the next decades radio surveys will be able to 
cover comparably wide areas with similar source number 
densities. In addition to that, the relative isolation of the 
few emission lines of astrophysical interest in the radio 
spectrum (e.g. the neutral hydrogen line at 1.4 GHz or 
molecular CO at 115 GHz) makes it possible to conduct 
intensity mapping observations, producing tomographic 
maps of the density fluctuations of these species. 

In this section we will describe two main cosmological 
probes of the low-redshift Universe in radio experiments: 
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intensity mapping (IM) of neutral hydrogen (HI) and 
continuum surveys. Both probes have complementary 
properties that make their cross-correlation extremely in¬ 
teresting: while intensity mapping is able to map the 
density field with remarkable radial (frequency) resolu¬ 
tion, but only on fairly big angular scales, continuum 
surveys lack any radial information, but are able to re¬ 
solve individual sources. On the other hand, while HI 
is known to be a very faithful tracer of the matter den¬ 
sity (i.e. has a low clustering bias), many radio sources, 
such as AGNs, are highly biased, and thus, as discussed 
in Section |HIB[ their combination could be a useful way 
to measure the level of primordial non-Gaussianity below 
the cosmic variance limit. 

For both probes we produce our forecasts assuming the 
current design for Phase 1 of the Square Kilometre Array 
(SKA) [3S1I35]. In both cases we assume that the exper¬ 
iment of choice will be SKAl-MID, an array of ~ 200 
single-pixel, 15m dishes to be installed in South Africa. 
It will cover the frequency range 350 —1760 MHz (z < 3) 
in two separate bands. Note that there are other inten¬ 
sity mapping experiments planned for the future, such as 
GHIME gZ] or BINGO [^. However, the SKA should 
provide a better coverage of the largest angular scales, 
which is the topic of this work. Additionally, the SKA 
is a multi-science facility that will perform both inten¬ 
sity mapping and continuum observations, and therefore 
cross correlating both datasets from the same experiment 
should be more straightforward. 


I. Intensity mapping 

While it is, in principle, possible to determine the 
radial position of radio sources through the redshift of 
the 21cm line, the intensity of this line from individ¬ 
ual sources is very small, and long integration times are 
needed in order to resolve it with sufficiently high signal- 
to-noise ratio. For a fixed amount of observation time, 
this limits the source number density, depth and area 
that can be achieved with such surveys. Intensity map¬ 
ping [4M5T] is a relatively novel technique that aims to 
circunvent this problem by measuring the combined emis¬ 
sion from all the sources in relatively wide pixels simulta¬ 
neously at different frequencies. Adding up the emission 
from hundreds or thousands of sources it is possible to 
detect the large-scale fluctuations in the HI density field 
in three dimensions up to relatively large depths. This is 
achieved at the cost of losing potentially valuable infor¬ 
mation on small angular scales, since individual sources 
are not resolved, although this is not necessarily an issue 
if we target large-scale observables, as is the case of this 
work. As was shown in m, however, two facts conspire 
against HI as a single-probe of non-Gaussianity and GR 
corrections: its relatively low clustering bias and the can¬ 
cellation of the linear perturbations on angular distances 
respectively. 

We model the clustering, magnification and evolu¬ 


tion biases for HI intensity mapping using the same ap¬ 
proach described in section 4.2 of m, based on using a 
scaling relation between HI and halo mass of the form 
M) (X M“, with a ~ 0.6. This allows us to esti¬ 
mate the HI density (and hence its evolution bias using 
equatio n [IT| ), as well as its clustering bias. As noted in 
section [II A[ the exact cancellation of the perturbations 
on transverse distances for intensity mapping implies an 
effective magnification bias shi = 2/5. 

As described in [S2], the most efficient use of SKAl- 
MID for intensity mapping will be in single-dish mode, 
and we base our noise model on this configuration. For 
an autocorrelation experiment, and assuming no correla¬ 
tions between the noise in different frequency channels, 
the noise power spectrum is constant on all scales: 

TVf-b ^ (31) 

where af is the noise variance per steradian in the i-th 
frequency channel, which can be related to the survey 
characteristics as 


2 ^ 47r /sky 

^tot A^dish 


(32) 


where ftot is the total observation time, /sky is the ob¬ 
served sky fraction, A^dish is the number of dishes in the 
experiment, Siy is the frequency bandwidth of the bin 
and Tsys = Tgky -I- Tlnst is the system temperature, which 
receives contributions both from the instrumental noise 
temperature and from the atmospheric and background 
noise (Tsky - 60K x (i// 300 MHzj-^-S). 

Regarding the radial binning, we follow the procedure 
used in m, subdividing the redshift range (0 < z < 3) 
into intervals of equal comoving width. We produced 
forecasts for gradually thinner bins until the uncertainties 
on our large-scale parameters, /nl and cgr converged. 
We thus obtain our final forecasts for 100 frequency bins 
with a comoving width Sx — 44Mpc//i. 

Since, in intensity mapping, sources are not detected 
individually, we must incorporate the effect of the tele¬ 
scope beam, and thus our model for the intensity map¬ 
ping power spectrum is 

Cf = BlBl + , (33) 


where is the harmonic transform of the instrumen¬ 
tal beam in the j—th frequency bin. For this work we 
have assumed that the beams are isotropic and Gaussian: 
Bg = exp{-£{£+l) 9%/2), with 9 3 ( 12 ) = c/(2.35 i^Ddish), 
where Ddish is the dish diameter. Note that when cross- 
correlating intensity mapping with any discrete tracer a, 
the effect of the beam is: 




(34) 


The experiment parameters used for SKAl-MID were: 
ri„st = 25K, /,ky = 0.75 (- 30000deg"), ftot = lO^h, 
-^dish ~ 200 Hnd- -/^dish ~ 15 m. 
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FIG. 6: Redshift distribution of radio sources for a flux cut 
S < 5/iJy. The distribution for the combined population 
is shown in dark grey, while the two subsamples considered 
here are plotted in blue (sample 1 - star-forming galaxies, 
starbursts and radio-quiet AGN) and green (sample 2 - FRI 
AGN). 


2. Continuum surveys 

The radio spectrum of most astrophysical sources is 
generally smooth and featurless, dominated by radio syn¬ 
chrotron emission with, perhaps, very few emission lines 
such as the aforementioned 21cm signal. It is thus possi¬ 
ble to integrate the flux from individual sources on rather 
wide frequency bands without losing too much informa¬ 
tion, while gaining access to much fainter sources than it 
would be possible to detect otherwise. Continuum sur¬ 
veys can therefore observe radio galaxies over incredibly 
large volumes with the caveat that all information regard¬ 
ing the radial distribution and clustering of these sources 
is completely unavailable. 

A number of papers have already shown the power 
of continuum surveys to constrain the level of primor¬ 
dial non-Gaussianity, especially making use of the multi¬ 
tracer technique m- Our aim here is to add on this 
result by considering the different cross-correlations with 
intensity mapping, as well as the potential of these cross 
correlations to detect relativistic effects. 

Our models for the redshift distribution and bias func¬ 
tions for the different radio sources considered here fol¬ 
low the method discussed in Appendix B of [TS], based 
on the measurements of the luminosity functions by [551 - 
Eg, used in [56] to produce semi-empirical simulations 
for radio surveys. We will consider here four different 
galaxy populations: normal star-forming galaxies, star- 
bursts, radio-quiet AGN and radio-loud AGN of the FR- 
I type. Owing to their very small number density, we 
omitted FRII galaxies from this analysis. Note that, as 
described in HU, separating the combined sample of ra¬ 


dio galaxies into these four different types could be a 
challenging task, especially given the enormous size of 
the sample that the SKA will be able to collect. Fur¬ 
thermore, although it should be possible to distinguish 
radio-loud AGN from other sources on a morphological 
basis, differentiating radio-quiet quasars would proba¬ 
bly require X-ray observations with a different experi¬ 
ment. Further separating normal star-forming galaxies 
from starbursts would require optical and redshift infor¬ 
mation, which would be difficult to obtain for the whole 
redshift range. For this reason we produce our forecasts 
for two cases: a pessimistic scenario with no source dif¬ 
ferentiation (i.e. one single tracer) and an optimistic one 
with two disjoint samples, given by radio-loud AGN on 
the one hand and a combination of star-forming galax¬ 
ies, starbursts and radio-quiet quasars on the other. We 
will call these two samples “sample 2” and “sample 1” 
respectively. 

As for intensity mapping, we consider here an SKAl- 
MID configuration, which should be able to detect radio 
sources to 0 ~ 5 over three quarters of the sky by inte¬ 
grating their continuum flux in the band 350-1050 MHz 
with an rms noise of 5'rms = I^JY) we will assume a 
5(7 detection threshold (i.e. ^cut = 5/iJy). The redshift 
distribution of the two populations considered here as 
well as the combined sample for this flux cut can be seen 
in Figure]^ and the corresponding clustering, magnifica¬ 
tion and evolution biases are shown in Figure together 
with those for HI intensity mapping. 

The noise power spectrum for continuum surveys is 
simply given by Poisson shot-noise where 

Nq is the angular number density of sources (per stera- 
dian) in the survey. We must note that the intensity 
mapping noise is not completely uncorrelated with that 
of continuum sources. At low redshifts most of the HI is 
located inside galaxies, and so many of the galaxies con¬ 
tributing to the HI line emission will also be detected in 
the continuum survey. However, the contribution of this 
Poisson term to the total noise power spectrum is very 
small compared to the instrumental noise levels |52j , and 
therefore we neglect these correlations. 


B. Forecasts 

We can expect the advantages of combining intensity 
mapping and continuum observations in a multi-tracer 
sense to be threefold: first, the large clustering bias of 
e.g. AGN combined with the very low b of neutral hy¬ 
drogen should significantly boost the amplitude of the 
/nl signal. Secondly, the absence of perturbations on 
transverse scales for intensity mapping should help to 
separate the GR effects and lensing magnification contri¬ 
butions in the cross-correlation. Finally, the main draw¬ 
back of continuum surveys is the absolute lack of radial 
information, and therefore, given the good redshift res¬ 
olution of intensity mapping, combining both datasets 
should significantly improve any constraints. 
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FIG. 7: Dependence of the constraints on /nl, cgr and ewL 
on a constant Gaussian prior on the clustering bias b{z). The 
constraints on ewL exhibit a much stronger dependence on 
Ab due to the fact that lensing magnification can only be 
constrained through the continuum survey. 


We have computed the forecasted errors on /nl, cgr 
and ewL for four different tracer combinations: 

• Case A: using only intensity mapping (IM only). 
These results were already presented in [iiliig. 

• Case B: combining intensity mapping and the joint 
continuum sample (IMxCont. (1 sample)). 

• Case C: combining intensity mapping and the two 
separated radio populations (IMxCont. (2 sam¬ 
ples)). 

• Case D: combining the two radio populations 
(Cont. only (2 samples)). 

In all cases we assumed complete overlap between the 
intensity mapping and continuum surveys (i.e. /sky = 
0.75). 

As in the case of LSST, we first look at the feasibility of 
detecting GR corrections by computing the uncertainty 
on ecR in the most optimistic case, with all parameters, 
except for /nl, fixed: 

o'A(eGR) = 2.75, crB(eGR) = 1-40, 
o'c(eGR) = 1-19, crD(eGR) = 10.4. 

Thus, even in the best-case scenario it would be impos¬ 
sible to observe the effect of the relativistic corrections 
above Icr. Following the argument laid out in section 
|III B| it is easy to understand why. First of all, the signal- 
to-noise level for both types of experiment is significantly 
lower than for photometric surveys. This is due to the 
smaller number density of radio sources above 5/iJy as 
well as the the high noise level expected for intensity 



FIG. 8: Joint constraints for /nl and ecR for the best three 
combinations of radio experiments: intensity mapping-only 
(purple), intensity mapping combined with a single continuum 
sample (grey) and intensity mapping combined with the two 
different continuum samples (green). Although the GR effects 
are not detectable, the improvement in the uncertainty on /nl 
due to the addition of continuum data is remarkable. 


mapping. Secondly, unlike in the case of red and blue 
galaxies, the magnification and evolution biases of the 
different tracers are not that different from one another. 

As was reported in m. the lack of radial information 
for continuum surveys makes the constraints on our large- 
scale observables very sensitive to the priors assumed on 
any nuisance or cosmological parameter. We can expect 
that the cross-correlation with intensity mapping should 
help in breaking some of these degeneracies, and there¬ 
fore we have studied them in detail. We obtain results 
that are mostly equivalent to those found for LSST: the 
constraints on cgr and /nl are insensitive to the prior 
on the clustering bias, but show a clear dependence on 
A/evo, although in this case they can be optimized with a 
looser prior of A/evo 1- The constraints on the ampli¬ 
tude of the lensing magnification term ewL, on the other 
hand, are significantly more sensitive to the priors on 
b{z) and s(z) than in the case of LSST. The main reason 
for this is that we can only constrain ewL through the 
continuum data, since the lensing term cancels exactly 
to first order for intensity mapping, and thus the lensing 
amplitude is more degenerate with other parameters. We 
show the dependence of ewL on Ab in Fig. 

As we did for LSST, we define an optimistic and pes¬ 
simistic set of priors, given by As = A/evo = 0.1 and 
As = A/evo = 1 respectively. In both cases we will as¬ 
sume that the clustering bias can be constrained to a 10% 
precision {Ab/b = 0.1) based on smaller-scale data [57j . 
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FIG. 9: Joint constraints on /nl, cgr and ewL for the combi¬ 
nation of both red and blue galaxies in LSST (black ellipses), 
for a combination of SKAl-MID intensity mapping with the 
full LSST sample (cyan) and for intensity mapping combined 
with both red and blue galaxies for (yellow). 


Our final constraints on /nl and ecR are shown jointly 
in Figure for optimistic priors and listed in Table [ll] 
Although the GR corrections remain unobservable for ra¬ 
dio experiments, the improvement of the constraints on 
/nl due to combining intensity mapping and continuum 
data is remarkable. 


VI. SYNERGIES 

In the last two sections we have described the poten¬ 
tial of two key next-generation experiments, SKA-MID 
and the LSST, to constrain cosmological observables on 
ultra-large scales, considering the potential multi-tracer 
analyses that could be performed using the capabilities 
of each experiment independently. We would now like to 
explore the possibility of combining the tracers probed by 
both surveys jointly, and the potential improvement that 
this synergy would represent over the results presented 
above. 

Both experiments will be located at similar latitudes 
and will therefore cover almost equivalent regions of the 
sky. There should also be a significant redshift over¬ 
lap between the LSST photometric survey and the SKA- 
MID HI intensity mapping experiment. The number of 
potential benefits of cross-correlating the datasets from 
both experiments is therefore immense, e.g.: the abil¬ 
ity of LSST to provide approximate redshifts for many 
of the sources in the SKA continuum survey, the ability 


of the SKA to provide more accurate redshift informa¬ 
tion through its HI survey to improve the LSST photo-z 
calibration, or the invaluable cross-checks that will be 
performed to mitigate the effects of survey-specific sys¬ 
tematic uncertainties (see [58] for a recent review of po¬ 
tential synergies). Since such potential benefits can not 
be ignored, this synergy will surely take place, and there¬ 
fore it makes sense to forecast for the potential impact 
on ultra large-scale cosmology. 

We will present our forecasts in two cases: 

• Combining SKA-IM and the LSST full sample, 
without differentiating between red and blue galax¬ 
ies. 

• Combining SKA-IM and the two red and blue LSST 
samples (i.e. three tracers in total). 

We do not consider cross-correlations between LSST and 
the SKA continuum survey for two reasons. Firstly, a 
large fraction of the radio galaxies detected by the SKA 
should be present in the LSST dataset too. Therefore 
these are not completely disjoint tracers, and evaluating 
their noise cross-power spectrum, crucial for the multi¬ 
tracer method, would be very uncertain. Secondly, the 
benefits of combining IM with a photometric survey are 
clearer: an improved coverage of the small radial scales 
due to the better spectral resolution of intensity mapping 
and an improved coverage of small transverse scales due 
to the better angular resolution of optical surveys. Fi¬ 
nally, as we discussed in section [V A H the low clustering 
bias, mild evolution bias and the absence of linear-order 
perturbations in the transverse distance hamper the abil¬ 
ity of intensity mapping to constrain /nl and cgr on its 
own. However, those same characteristics make it a very 
useful probe to combine with other tracers in the multi¬ 
tracer scheme. 

Although both experiments should cover almost the 
same area of the sky, we can expects a certain loss of over¬ 
lap caused by regions of the sky dominated by different 
systematics (e.g. dust extinction for LSST and galactic 
synchrotron for SKA-IM). The forecasts presented here 
conservatively assume that both experiments will yield 
useable observations in a patch of area /sky = 0.4. Note 
that we treated this value as the sky fraction in which 
both experiments can be cross-correlated, however we 
still use the full SKA area (/sky = 0.75) when computing 
the noise power spectrum for intensity mapping. 

The full LSST sample was modelled by combining the 
red and blue luminosity functions, and conservatively as¬ 
suming a photo-z uncertainty erg = 0.05. Following the 
prescription outlined in Section [IV A[ this entails divid¬ 
ing the sample into the same 9 redshift bins used for the 
blue population. 

The joint constraints on /nl, cwl and cqr for the two 
cases above, together with the constraints found for LSST 
alone are shown in Figure where an optimistic prior 
As = A/evo = 0.1 was used. The quantitative results are 
listed in Table [III for both pessimistic and optimistic pri¬ 
ors. Several conclusions can be drawn from this exercise: 
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• In spite of the low clustering bias of neutral hydro¬ 
gen, the errors on /nl improve immensely by com¬ 
bining LSST and intensity mapping (by a factor of 
~ 3 — 4 compared to the LSST-only case). 

• Although we observed that the large evolution bias 
of red galaxies was crucial in boosting the detection 
level of GR corrections when jointly analyzing the 
red and blue samples for LSST, we now see that 
the cross-correlation of HI intensity mapping and 
the full LSST sample, which does not display such a 
large evolution bias, could still yield a measurement 
of cgr above the 5a level. In the most optimistic 
case (optimistic priors and a joint analysis of the 
three samples) this increases to ^ 20a. 

• Due to the cancellation of the lensing term for in¬ 
tensity mapping the improvement in the detection 
of the large-scale lensing magnification is only mild. 


VII. DISCUSSION 

In this paper we have shown how more than one tracer 
of the cosmic density field can be used to extract pre¬ 
cise measurements of the ultra large-scale features of the 
universe. We have shown that a judicious combination 
of photometric redshift surveys with cosmological radio 
surveys will allow us to obtain precise constraints of pri¬ 
mordial non-gaussianity, of weak lensing on large scales 
and most significantly, of general-relativistic effects which 
depend on various combinations of the gravitational po¬ 
tential. Our final results are summarized in Table El in 
terms of the marginalized uncertainties on /nl and the 
amplitudes of the GR corrections (cgr) and lensing mag¬ 
nification (cwl)- 

More specifically, we have shown that the best way to 
optimize the detection of GR effects is to use a combina¬ 
tion of low-noise tracers (i.e. high number density) with 
very different evolution and magnification biases. We 
find that ongoing and future photometric redshift sur¬ 
veys will be able to satisfy both conditions, given the 
large expected number density of observed sources and 
the possibility of separating the galaxy sample into two 
populations, blue and red galaxies, with different evolu¬ 
tion of their luminosity functions. Indeed we have seen 
how the rapid decay in the number density of red galaxies 
above redshift ~ 1 translates itself into a large /evo that 
boosts the signal of the GR corrections and could enable 
LSST to detect them with a significance of ~ 5—10 a. Us¬ 
ing the same method we have also shown that DES might 
be able to make a lower significance detection (~ 3 a) in 
a shorter term. 

We have also explored the possibility of combining two 
of the main cosmological radio probes, HI intensity map¬ 
ping and continuum surveys, planned for the Square Kilo¬ 
metre Array. We have seen that in this case, the higher 
noise levels make a detection of the GR corrections unvi- 
able, in spite of the larger volume covered in comparison 


with LSST. However, the contrast between the low clus¬ 
tering bias of neutral hydrogen and its large value for cer¬ 
tain radio sources makes the combined analysis of both 
probes very beneficial to reduce the uncertainties on the 
level of primordial non-Gaussianity. We have seen that, 
in the worst-case scenario (1 single continuum sample and 
pessimistic priors), the error on /nl could improve by a 
factor of 3 with respect to the constraints from intensity 
mapping alone. 

Finally, since LSST and SKA will carry out their obser¬ 
vations over a very similar region of the sky, the benefits 
of cross-correlating both experiments will be immense for 
many science cases, and here we have studied the poten¬ 
tial impact of this synergy in the study of cosmology on 
ultra-large scales. We have shown that it would be pos¬ 
sible to improve the constraints on /nl by a factor of 
at least 2 — 3 with respect to the individual constraints 
from each experiment. The detection level of the GR 
effects would also improve significantly and, more inter¬ 
estingly, the combination of the full LSST sample and 
SKAl intensity mapping should also be able to detect 
these effects above 5a regardless of the large evolution 
bias of red galaxies. 

We have shown that the detectability of these new ef¬ 
fects is an exciting prospect, yet many important chal¬ 
lenges must be overcome before such measurements can 
be carried out. Firstly, the estimates of clustering statis¬ 
tics on ultra-large scales in most experiments are rid¬ 
dled with systematic uncertainties. Intensity mapping 
is a perfect example of this: galactic and extragalactic 
foregrounds have radio intensities that are several or¬ 
ders of magnitude larger than the cosmological HI sig¬ 
nal. Although both components can in principle be sep¬ 
arated based on their very different spectral behaviour, 
any foreground residual will necessarily contaminate the 
signal on large scales, both angular (due to the large- 
scale features of galactic foregrounds) and radial (due 
to their smooth frequency dependence). As shown in 
m, imperfect foreground cleaning can lead to the cross¬ 
correlations between well separated frequency bins being 
unreliable, and, as was reported in m, excluding these 
cross-correlations from the analysis can have a very detri¬ 
mental effect on the uncertainty on large-scale parame¬ 
ters such as /nl- Other instrumental issues, such as the 
effect of correlated 1// noise in the data, a bad charac¬ 
terization of the instrument’s beam or the contamination 
from polarized foregrounds due to polarization leakage 
can introduce even more troublesome systematics in the 
data that must be correctly minimized in order to make 
these measurements possible. 

Optical surveys are also subject to a sizeable collection 
of systematics. Effects such as galactic dust extinction 
and airmass can modify the observed magnitude and gen¬ 
erate a non-uniform magnitude limited sample which can 
modify the observed clustering statistics on large scales. 
Star contamination, i.e. the accidental inclusion of stars 
in the galaxy sample, and star obscuration, as reported 
in IE], can also have a significant impact on the clus- 
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Experiment type 

Tracers 

<^(/nl) 

cr(eGR) 

cr(ewL) 

Photometric survey 

LSST, red-only 

4.53 (4.54) 

1.65 (1.70) 

0.18 (0.67) 

(LSST) 

LSST, blue-only 

1.71 (1.72) 

3.45 (3.48) 

0.12 (0.61) 


LSST, red x blue 

1.62 (1.63) 

0.10 (0.17) 

0.10 (0.36) 


DES, red x blue 

7.18 (7.20) 

0.29 (0.32) 

0.06 (0.06) (*) 

Radio 

IM-only 

3.00 (3.01) 

2.71 (2.75) 

— 

(SKAl-MID) 

IMxCont., 1 sample 

0.86 (0.89) 

1.47 (1.95) 

1.16 (3.80) 


IMxCont., 2 samples 

0.69 (0.71) 

1.23 (1.43) 

0.33 (2.21) 


Continuum-only, 2 samples 

1.91 (1.97) 

23.9 (31.1) 

0.56 (2.85) 

Synergy 

IMxall 

0.41 (0.41) 

0.15 (0.39) 

0.09 (0.32) 

(SKAl-MID X LSST) 

IMx red X blue 

0.40 (0.40) 

0.05 (0.12) 

0.07 (0.19) 


TABLE II: Forecasted constraints on /nl, egr, and ewL for the different experiments explored in this work. In each case, the 
constraints assuming optimistic and pessimistic priors on the magnification and evolution biases are shown with and without 
parentheses respectively. (*)Note that the uncertainty on ewL is shown to be smaller for DES than LSST. This does not mean 
that DES would be better at measuring lensing, only that the value of the magnification bias at the magnitude limit of DES 
(~ 24) enhances the lensing magnification effect. LSST would be able to greatly improve on this result by adopting a different 
(lower) flux cut. This was already noted in [59] . 



FIG. 10: Dependence of the constraints on cgr (solid) and 
/nl (dashed) on the largest angular scale included in the anal¬ 
ysis (given by the minimum multipole ^min), for three different 
tracer combinations: LSST red and blue galaxies (black), in¬ 
tensity mapping and our 2 continuum samples (green) and IM 
combined with both LSST samples (yellow). The constraints 
degrade significantly after £min ~ 10. 


tering observed on large angular distances, since both 
effects grow towards the galactic plane. Another specific 
source of systematic uncertainty for photometric surveys 
is the incorrect characterization of the window functions 
corresponding to each redshift bin. This can be caused 
by uncertainties in the galaxy redshift distribution, due 
to the absence of precise redshifts, and to inaccuracies 
in the characterization of the photo-z distribution as a 
function of redshift. 

Most of these systematic effects, unless corrected for, 


will limit our ability to use the observations on the largest 
angular scales. We have studied the dependence of our 
forecasts on the largest scale included in the analysis (i.e. 
the minimum angular multipole fmin) in Figure for 
our best three tracer combinations (LSST red and blue, 
IM combined with the two continuum samples and IM 
combined with the two LSST samples). The constraints 
on /nl and ecR degrade significantly above £min ~ 10. 

In addition to these systematic effects, we have also 
identified other challenges that must be tackled before 
this detection can be achieved. We have found that hav¬ 
ing at least two tracers with very different magnifica¬ 
tion and/or evolution biases is key to boosting the signal 
of GR effects, and we have illustrated this point with 
the case of red and blue galaxies in photometric surveys. 
However, we have also seen that, even if such tracers can 
be found, a sufficiently precise knowledge of their joint 
magnitude-redshift distribution is necessary to statisti¬ 
cally detect the GR effects, which would otherwise be 
swamped by the uncertainties on s(z) and /evo- In the 
case of LSST we have shown that priors of the order of 
As ~ A/evo — 0.1 would be needed to optimize this 
measurement. To put these uncertainties into perspec¬ 
tive, we should note that the uncertainty in the value of 
/evo propagated from the errors in the estimate of the 
luminosity function for red galaxies made by m, based 
on data from the DEEP2 and GOMBO-17 surveys, is 
of the order of A/evo ~ 2 at high redshifts {z > 0.8), 
where the GR effects are maximally boosted. Assuming 
that this uncertainty decreases with the square-root of 
the observed number of galaxies, spectroscopic observa¬ 
tions of a representative sample of ~ 1.5 million galaxies 
with redshifts above ^ 0.8 would be needed to achieve 
the optimal sensitivity. This would correspond to a frac¬ 
tion of ~ 2% of the total red sample expected for LSST in 
the same redshift range. The Dark Energy Spectroscopic 
Instrument (DESI) |Slj, for instance, should be able to 
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measure spectra for ~ 4 million luminous red galaxies in 
the redshift range 0.4 < z < 1.0, although these observa¬ 
tions would be carried out in the northern hemisphere. 

Finally, it is worth noting that the detection levels and 
uncertainties reported here depend critically on the fidu¬ 
cial models used for the bias functions of the tracers un¬ 
der study. Although we have tried to use reliable and 
conservative assumptions for these models, these quan¬ 
tities are still uncertain, and therefore the quantitative 
results could change as those uncertainties are reduced. 
On the other hand, this also implies that these results 
could actually be improved if one could find more op¬ 
timal tracer classifications (e.g. beyond the traditional 
red-blue split) enhancing the amplitude of the relativis¬ 
tic corrections. 

If we can tackle all these challenges, then, by measur¬ 
ing ultra large-scale modes in the density field, we will 
be able to break out of what has been the conventional 
arena of large scale structure studies - the quasi-static 
regime - on substantially sub-horizon scales. Until now, 
most surveys of large-scale structure (such as redshift or 
weak lensing surveys) have probed what can be thought 
of as the Newtonian regime of cosmological perturbation 
theory, where the time dependence of the gravitational 
potentials can be ignored [63] . Now, by probing scales 
which span the cosmological horizon, this time depen¬ 
dence comes into play, giving us a completely new win¬ 
dow on the evolution of cosmic perturbations. 

Accessing ultra-large scales can be of particular signifi¬ 
cance in constraining modifications to General Relativity 
[26| . It has been shown that in the quasi-static regime, 
cosmological data will only constrain a limited subspace 
of possible modifications to gravity [64ll66j . To extend 
the scope of the constraints, it is essential to break any 
degeneracies that arise on small scales by using the time 
variation of the gravitational potentials on large scales. 
As shown by |67||68] one will expect specific signatures (of 
order of a few percent) which, with the methods proposed 
in this paper, might be observable with future surveys. 

Given that we are targeting ultra-large scales, one 
should bear in mind the mild evidence for deviations from 
the standard, AGDM, concordance model on large scales. 
The recent Planck measurements of the the GMB on large 
scales has found a number of anomalies - suppression 
of power on very large angles, correlations between the 
large angle multipoles and an anomalously cold spot - 
which have yet to be understood |69j . Gosmic variance 
severely restricts any attempts at better characterizing 
these anomalies with the GMB alone. One would hope, 
however, that with the possibility of better constrain¬ 
ing large scale modes with multiple tracers, one could 
substantially improve our characterisation of these large 
scale anomalies with complementary constraints. 


We have focused on extracting information at the level 
of auto and cross-correlations of the various tracers. An 
interesting, uncharted, avenue would be to attempt to 
use these techniques to extract maps of each of the effects 
that come into play. In other words, and much like what 
is done in the case of the GMB HQ], one could imagine 
constructing linear combinations of the various observ¬ 
ables (which are three dimensional maps of a particular 
tracer) to extract, for example, maps of the gravitational 
potential, their derivatives, large scale flows, etc. This 
is a somewhat ambitious goal but if achievable would 
open up a completely new window into the large scale 
properties of space-time and matter; a direct compari¬ 
son between the different maps, and a comparison be¬ 
tween these maps and other cosmic observables might be 
used to uncover different cosmological properties which 
are usually buried in the standard correlation function 
(or power spectrum) analysis. 

Finally, we have focused on a particular choice of 
current and future surveys - photometric redshift sur¬ 
veys and cosmological radio surveys - but this does not 
cover all possible observational approaches. Spectro¬ 
scopic galaxy redshift surveys, such as Euclid lu and 
WFIRST [72], will also attempt to constrain the effects 
we discuss in this paper, albeit on smaller scales and 
with a lower source density. A more speculative, ambi¬ 
tious and, in our mind, exciting approach using an all sky, 
wide band survey, SPHEREX, would have all the right 
characteristics for the methods proposed here. As ar¬ 
gued in the number density and bias of these sources 
should be sufficiently large to aid, and complement, the 
constraints that will be found with the combination of 
probes found here. In addition to these probes of galaxy 
clustering, it would also be worth exploring the poten¬ 
tial benefits of combining them with other cosmological 
observables, such as weak lensing, GMB lensing or even 
the GMB itself. These combinations could benefit from 
the particular properties of these tracers (e.g. zero bias 
for weak lensing and minimal clustering variance for the 
GMB), significantly improving the constraints reported 
here. We are clearly entering an era in which precision 
measurement of ultra large-scale modes can become a 
core science goal of observational cosmology. 
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Appendix A: The linear order components of the 
fluctuation in source number counts 

The complete expressions, to linear order, for the 
transfer function of the perturbation in the source num¬ 
ber counts used in Equation]^ is given by [20] : 


A^’\k) = J dr]bW^6M,syn{k,r])ji{kx{r])), Af^^’\k) = JdiqiaH) e{k,'q) jg {kx{r])), (Al) 

A^’\k) = £{£ + 1) J dr]Wl^{r]){(j)-\-'ip){k,r))je{kxiv)), Aj'^'\k) = J drj (/evo - 3) aHWi{r]) jeikxiv)), 

(A2) 

= J drj (^1-^ -k + 5s - /evo^ W^{r]) feikxiv)), (AS) 

A^^'\k) = J dr] (^2-\--\-5s-fe^o'^W,{r])ip{k,r])ji{kxiv)), (A4) 

A^^’\k) = J df]{5s-2)W^{r])(j){k,7])ji{kx{'n)), A^^’"{k) = J dr] {aH)-^W,{r]) (j)'{k,r]) ji{kx{v)), (A5) 

A^^’\k) = J dr]Wf’^{r]){(t) + ip){k,r])je{kx{v)), Af^'^k) = J dr]Wl^'^{r]){(j) + i;y{k,r])je{kx{v)), (A6) 


where ji{x) is the spherical Bessel function of order £ and where we have defined the window functions 


Wi{r]{z)) = Wi{z) 
ru 


( ^ 

V dz 


wriv) ^ 


9 — 

dv'Wyi]') 


wtiv) = I dr]'Wi{r]') 


1,2- 5s(r]') x(v) - X(v') 


x(v)x(v') 


/o 


X 




+ 5S - /e- 


The quantities SM,syn, d, (j) and '0 in the equations above 
are respectively: the transfer function for the matter den¬ 
sity perturbations in the comoving synchronous gauge, 
the transfer function for the peculiar velocity divergence 
in the conformal Newtonian gauge and transfer functions 


for the two metric potentials in the same gauge, defined 
by the line element: 

ds^ = —a^(r]) [(1 -k 2'ij})drf — (1 — 2(j>)Sijdx'^dx^ . (A7) 














